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At very high energies, the relevant component of the hadron wavefunction can be 
described as a Color Glass Condensate, i.e., a state of high density gluonic matter 
whose distribution is random, but frozen over the relevant time scales. The weight 
function for this distribution obeys a renormalization group equation in the form 
of a functional Fokker-Planck equation. Its solution leads to an effective theory 
which predicts gluon saturation at sufficiently high energy, or small Bjorken's x. 



1 Introduction 

Hadronic scattering at high energy, or small Bjorken's a;, uncovers a novel 
regime of QCD where the coupling is small (a^ ^ 1) but the parton densities 
are so large that cmiyentioiLal perturbation theory breaks down, via strong 
non-linear effects Ero'B (see CrO for recent reviews and more references). Re- 
markable progress in recent years has resulted from the observation that this 
high density regime can be studied via semiclaasicaLmethodsQ. In what fol- 
lows, I shall describe a classical effective theoryQ^Q'LrO which is well suited to 
study the non-linear phenomena at small x and has a transparent physical 
interpretation: It portrays the gluon component of the hadron wavefunction 
(the relevant component at small x) as a Color Glass Condensate (CGC). 
This is a multiparticle quantum state with high occupation numbers, but to 
the accuracy of interest it can be represented as a stochastic color field with 
a probability law determined by a functional Fokker-jElanck equation. 

The latter is a renormalization group equationErQ (RGE) which shows 
how to construct the effective theory by integrating out quantum fluctuations 
in the background of a strong color field (the CGC). The non- linear effects 
included in the RGE via this background field describe interactions among the 
gluons produced in the quantum evolution towards small x. This leads to a 
non-linear generalization of the BFKL equation which, remarkably, predictffl 
the saturation of the gluon distribution at small x : Unlike the BFKL result 
for the gluon distribution, which increases exponentially with the rapidity 
T = ln(l/a;) cx Ins and thus violates the Froissart unitarity bound a < In^ s, 
our corresponding prediction grows only linearly with t, thus being consistent 
with the unitarity bounds. 



2 The effective theory for the CGC 



The effective theory applies to giuon correlations in the hadron wavcfunction 
as measured in deep inelastic scattering at small Bjorken's x. It is formulated 
in the hadron infinite momentum frame, where small x corresponds to soft 
longitudinal momenta^ fc+ = xP^, with P+ the hadron momentum, and 
a; <C 1. The main observation is that the "fast" partons (i.e., the excitations 
with p"*" ^ fc"*") can be replaced, as far as their effects on the soft correlation 
functions are concerned, by a classical random color source p"^ (x) , whose gross 
properties are determined by the soft-fast separation of scales: 

The fast partons appear to the soft gluons as a color charge distribution 
p°'{x) which is static (i.e., independent of a;^), localized near the light-cone 
(within a small distance Aa;~ ~ *C l/fc"''), and random (since this is the 
instantaneous color charge in the hadron "seen" by the soft gluons at the arbi- 
trary time of their emission) , with gauge-invariant probability density Wr [p] . 
(We use the rapidity r = ln(P+/A:"'") = ln(l/a;) to indicate the dependence of 
the weight function upon the soft scale .) 

Gluon correlations at the scale k'^ = xP^ are obtained as {x = {x^ , x±)): 

{Al, ix+ , x;)Al {x+,y)---)r^ JvpWr [p] Al, ix)Al (y) • ■ • , (1) 

where Al^ = A^ [p] is the solution to the classical Yang-Mills equations 

(D,P''^)4a;) = S^'+paix) (2) 

in the light-cone gauge A'^ = 0, which is the gauge which allnws for the most 
direct contact with the gauge-invariant physical quantities □'Q. For instance, 
the gluon distribution function is obtained as 

xGix,Q')^^J |^e(Q2-fci)(|^+^(fc)p)^. (3) 

where = d^A^^ is the electric field associated to the classical solution 
Al[p], and k = {k+,kj_) with fc+ = xP+ = P+e"^. 

Eqs. (|l|) and (||) are those for a 17/055 (here, a color glass): There is a 
stochastic color charge, that is averaged over. This is frozen over short time 
scales, of the order of the lifetime Aa:+ ~ l/k~ oc fc+ of the soft gluons, but 
it changes randomly over the larger time scale l/p^ ^ l/fc , which is the 
natural time scale for the dynamics of the fast partons. This is like a glass 
which is disordered and is a liquid on long time scales but seems to be a solid 

"^I use light-cone vector notations, e.g., fc'' = (k'^ , k~ , k j_) , with A:+ = (fc" + k^)/\/2, 
k- = (fc" - fc^)/V2, fcx = {k^,k^), and fc ■ x = k~x+ + k+x~ -k^-x^. 



on short time scales. We shall see later that, at saturation, the gluons are very 
densely packed, with a phase-space density ^ which is the maximum 

density allowed by their mutual interactions. They form a Bose condensate. 

In the effective theory, this Color Glass Condensate is described by strong 
classical color fields Al^ ~ The classical dynamics is then fully non-linear, 
so one needs the exact solution to eq. (||), which readsEl 

A'ix) = -U{x)d'U^ix), (4) 

J — oo 

-Via(f) = p{x). (6) 

Since p{x), and therefore a{x), are localized near a;~ = 0, so is the electric 
field, which appears effectively as a (5-function to any probe with low longitu- 
dinal resolution (i.e., with momenta < fc+) : 

T+\x) ^ 6{x-)-Vid'V''), V''ix±_) = Pcxp\^ig f dx- a"" (x' , x ^ jT^Y (7) 



3 The Renormalization Group Equation 

The weight function WV[p] for the effective theory at the scale A;+ = xP~^ 
is obtained by integrating out the quantum fluctuations with p+ > fc"*" in 
layers of p^, to "leading logarithmic accuracy" — i.e., by retaining only the 
terms enhanced by the large logarithm ln(l/a;) — , bU|L|±o all orders in the 
background fields .4^ generated at the previous steps H'B. The background 
fields are an essential ingredient: their interactions with the quantum gkions 
correspond, in the more conventional picture of the "parton cascades" El, to 
rescatterings among the soft gluons radiated in the quantum evolution towards 
small x (e.g., recombinations of gluons from different cascades). 

The quantum evolution can be formulated as a renormalization group 
equation (RGE) for the flow of WV[p] with t — ln(l/a::). To motivate the 
structure of this equation, I succintly describe one step in this RG procedure. 
Assume that we know the effective theory at some initial scale A"*" — as 
specified by Wa[p] = W^t[p], with r — ln(P+/A+) — , and we are interested 
in correlations at the softer scale ~ 6A+ with 6 <C 1 and asln{l/b) < 1. 
Our purpose is to construct the new weight function WhA[/o] = WV-i-At[p] 
(At = ln(l/5)), which would determine the gluon correlations at this softer 
scale. As compared to W\ [p] , the new weight function must include also the 



quantum effects induced by the "semi-fast" gluons with p+ in the strip 

5A+ < \p+\ < A+. (8) 

To characterize these effects, consider the change in the 2-point function 
{Al^{x)Al^{y))\ with decreasing A"*". At the original scale A+, we can use 
the effective theory to write {A'{x)A'{y))A = {A'{x)A'{y))wA: with A'jp] 
the classical solution in eq. and the average over p computed as in eq. @). 
At the new scale bA~^, AJ, = Al^ + SAl^, with 5Al^{x) representing the quantum 
fluctuations with momenta < A+. Thus, 

{A\x)A\y))^^ = ({{A' + SA%{A' + 6A%) ) 

^{A\x)A\y))^^^, (9) 

where the internal brackets (■ • ■)p in the first line stand for the quantum 
average over semi-fast fluctuations (cf. eq. (||)) at fixed p, while the exter- 
nal brackets (■ • ■)vKa denote the classical average over p with weight function 
Wa[p]- By definition, the result of this double averaging (i.e., of the classical 
plus quantum calculation) in the original effective theory at the scale A"*" must 
be the same as the result of a purely classical calculation in the new effective 
theory at the scale foA"*" . This is the content of the second line in eq. (||) . 

The quantum corrections in the first line of eq. must be computed 
to lowest order in Q!sln(f/6), but to all orders in the background fields 
This is essentially an one-loop calculation, but with the exact background 
field propagator of the semi-fast gluons. By matching its result with the clas- 
sical calculation in the second line in eq. (||), one can deduceCl the functional 
change ISW = Wr+Ar — Wr necessary to absorb these new correlations. Since 
AW (X At, the evolution is formulated as a (functional) RGB, which is most 
conveniently written for WV[a] = Wr[p — — V^^a], cf. eq. (^). It readsQ|^ : 

(10) 

where riab{x±, y±) it itself a non- linear functional of the color field Uaix), via 
the Wilson lines V and defined in eq. (|^) : 

(11) 

''In Ref.i, Weigert has shown that Balitsky's equation can be summarized into a func- 
tional equation equivalent to (|lo|. More recently, Sr-simplified derivation of eq. ( p^ ) has 
been given by Mueller, from the dipole point of viewEJ- 



The functional derivatives in eq. ( pj]| ) are taken with respect to a'!^{x±) = 
a°'{x~ = x~,x±), where x~ = 1/A+ ^ e"^. This is so since the quantum 
corrections are located in the strip 1/A+ ^ x^ ^ 1/6A+, i.e., on top of the 
original field which has support at < ^ 1/A+. Thus, by integrating 
out the quantum modes in layers of p"*" , one constructs the classical field (or 
source) in layers of x~ , with a one-to-one correspondence between and x~ 
which reflects the uncertainty principle A.T^Ap+ ~ 1. 



4 General properties and consequences 

Eq. ( p^ ) has the structure of a diffusion equation: It is a second-order (func- 
tional) differential equation whose r.h.s. is a total derivative, as necessary to 
conserve the total probability. It describes quantum evolution as the diffu- 
sion (with "time" r) of the probability density Wr [a] in the functional space 
spanned by aa{x~ ,x±). 

For comparison, consider the usual Fokker-Planck equation describing 
Brownian motion in flat space {P{x, t) = the probability density to find the 
particle at point x at time t knowing that it was at a; = at time 0) : 

dP{x,t) ^ d 



dt 



where D is the diffusion constant and F' — —dV/dx^ is an external force. If 
= 0, the corresponding solution: 

(iJo^M-i^} 

goes smoothly to zero at any x when t ^ oo (runaway solution). But for 
1^ 7^ 0, there exists a non-trivial stationary solution Pq{x) ^ exp[— y(a;)/Z)], 
i.e., an equilibrium distribution that is asymptotically reached by the system. 

There is manifestly no force term in the RGE (p^. This jpgoperty relies 
on subtle compensations between real and virtual corrections EI'LJ, and entails 
a purely diffusive behaviour, as in eq. (^). Thus, there is no fixed point, or 
stationary distribution, for the flow described by eq. (pO|). 

From the functional RGE (p^), one can derive ordinary evolution equa- 
tions for all the observables which can be computed as an average over a . 
Particularly interesting quantities are the gauge-invariant products of Wilson 
lines {iT:{V,}Vy)tT{y^Vu) ■ ■ •)r, for pszhich eq. ( |lO|) predicts the same evolution 
equations as obtained by Balitskyll3 via operator product expansion. In par- 



ticular, the 2-point function {tr{VjVy))r satisfies 



X {NMVlVy) - tr(FjT4)tr(V;V,))^ (14) 

which in the large Nc limit becomes a closed equation, due to KovchegovEl, 
that can be recognized as a non- linear generalization of BFKL. Its solution 
determines the cross section ddipoie i'^mf^ ) for the scattering of a "color dipole" 



of size r^=x^— y± off the hadron cIeJI. Recent progresM^with this equation 
exhibits color transparency {(Jdipoie oc r^Qlir)) at small distances rj^ <C 
1/Qs(t), and saturation {adipoie ~ const) at large distances r±_ ^ 1/Qs(t), 
in qualitative jagreement with a phenomenological model by Golec-Biernat 
and Wiisthoff Il3 which successfully describes the HERA data. In the above 
estimates, Qs{t) is the saturation scale to be discussed in the next section. 



5 Gluon Saturation and Unitarity 



Except in the large Nc limit, Balitsky's equations do not close individually, 
but rather form an infinite hierarchy of coupled equations. It is then mope 
convenient to study directly the functional RGE (^^. As shown in RefJa, 
approximate solutions to this equation can be obtained in limiting kinemati- 
cal regimes, by combining mean field approximations (MFA) and kinematical 
simplifications. Moreover, an exact solution_in the form of a path-integral 
in 2+1 dimensions, can be also writtejijdowntfl, which may be used in lattice 
calculations. (See also the second RefJlJ for a different path-integral formula- 
tion.) Here, I shall briefly describe the analytic solutions, with emphasis on 
the phenomenon of gluon saturatiorjj. 

The MFA consists in replacing riab{x±,y±) in the kernel of the RGE ( |l0| ) 
by its expectation value {'riab{x±,y±))T = <5°^7r (a;_L, j/±). Then, the equation 
can be solved exactly, and the solution Wr — a Gaussian in a (or p) — can 
then be used to self-consistently determine the approximate kernel 7r(a;±, y±)- 

The kinematical approximations rely on the fact that there is an intrinsic 
transverse scale in the problem, the saturation scale Qs (t) , which is the inverse 
correlation length for Wilson lines: 

\0, for -y^l > 1/Q^(t) ^ 

To get an estimate for Qs{t), one starts the analysis at sufficiently short 
distances \x± — y±\ <C 1/Qs{t), where the fields are weak and their dynamics 



is linear {Vj « 1 + iga{x±)), and study the onset of non-linearities with 
increasing \x± — yj_\. In this regime, one obtains (cf. eq. @): 

(V:Vy)r ^ exp{ - ^^xG{x,l/rl)}, (16) 

(cis = asN/ir, r± ~ x± — yj_, and R the hadron radius) with xG{x,Q^) 
satisfying the standard evolution equation in the double log approximation: 

^ ^2 Q') = asxG{x, Q2). (17) 

OT a log y ^ 

(Here, this equation arises via the self-consistency requirement.) Non-linear 
effects, which are negligible at sufhciently small r^, become important for 
rj^ ~ 1/Qs{t) where the exponent in eq. ( |l^ ) becomes of order one. Thus, 

R^QHt) ^ asxG{x,Ql{T)), (18) 

which together with eq. ( [l7| ) implies that Qs{t) increases exponentially: 

QsW = g^(ro)e'''='"~"°^ (c = 4inDLA). (19) 
The weight fmiction at high- momenta kj_ 3> Qs{t) is the Gaussian: 

[,] ^ expf - 1 / P'ik^)P^i~k.) \ ^ (20) 



(2^)3 ^,(fc^) 

where ^r{k±) oc (9/(?log fc^)a::G(a;, A;^) is only slowly varying with 
Eq. ( pO| ) is the McLerran-Venugopalan model of independent colour charge 
amended by the standard, linear, quantum evolution, cf. eq. (p7|). 

The new physics, intrinsically non-linear, shows up at large distances 
(r_L > 1/Qs(t)), or low momenta {k± < Qs^t), with ki_ > Aqcd though), 
where color fields are strong, a° ~ 1/g, and Wilson lines are decorrelated, 
Vj Vy ft ! 0. Then one can neglect the bilinears involving Wilson lines in 
eq. (|11[), so that the corresponding weight function is very simple: 



Wi' 



-\p] - expl - i / ^ p°(fcJp°(-fcJ l (21) 

This describes a 2-dimensional Coulomb gas, with "dielectric constant" 

er{k^) = r - f(fci) = ^ \n^P-, (22) 

which is the rapidity window for quantum evolution in the saturation regime. 
Indeed, for given fcx ^ Qs (''"), only those partons are saturated whose rapidi- 
ties y are large enough for Qs{y) > k±. This requires f{k±) < y < t, with 
f{k±) the rapidity where Qs(T(k±)'^ = k±. 



Saturation is an immediate consequence of eq. (pT|), which impHes the 
following gluon density per unit transverse phase space at k±_ < Qs{t) '■ 

At fixed r, the gluon density is almost constant]^ for fc^ < Qs{t)^ and of order 
0{\/as)- This is the Color Glass Condensate. When r ~ Ins increases, the 
gluon density at < Qs{t) increases only linearly with t (i.e., logarithmi- 
cally with the energy s), but the saturation scale Qs{t) itself increases as a 
power of s. That is, with increasing energy, the new partons are predom- 
inantly produced at large transverse momenta ^ Qs{t), so they cannot be 
discriminated by an external probe with resolution < Q'^{t). Thus, al- 
though the total number of gluons keeps growing when the energy increases, 
there is no contradiction with unitarity since, beyond some energy, the par- 
tons produced by quantum evolution are too tiny to contribute to the cross 
section at fixed Q^. Saturation is a natural mechanism to restore unitarity. 
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